EPJ manuscript No. 

(will be inserted by the editor) 



Nuclear mass form factors from coherent photoproduction of 



7T° mesons 



in 
o 
o 

(N 
Oh 

c/3 



> 

o 
o 

o 
in 
o 

=5 



X 



B. Krusche 

Department of Physics and Astronomy, University of Basel, Ch-4056 Basel, Switzerland 
Received: date / Revised version: date 

Abstract. Data for coherent photoproduction of it mesons from nuclei ( 12 C, 40 Ca, 93 Nb, na *Pb), recently 
measured with the TAPS detector at the Mainz MAMI accelerator, have been analyzed in view of the 
mass form factors of the nuclei. The form factors have been extracted in plane wave approximation of 
the A('y,ir°)A reaction and corrected for final state interaction effects with the help of distorted wave 
impulse approximations. Nuclear mass rms-radii have been calculated from the slope of the form factors 
for q 2 — » 0. Furthermore, the Helm model (hard sphere form factor folded with Gaussian) was used to 
extract diffraction radii from the zeroes of the form factor and skin thicknesses from the position and height 
of its first maximum. The diffraction radii from the Helm model agree with the corresponding charge radii 
obtained from electron scattering experiments within their uncertainties of a few per cent. The rms-radii 
from the slope of the form factors are systematically lower by up to 5% for PWIA and up to 10% for 
DWIA. Also the skin thicknesses extracted from the Helm model are systematically smaller than their 
charge counter parts. 

PACS. 13.60.Le meson production - 25.20.Lj photoproduction reactions - 21.10.Gv mass distributions 



1 Introduction 

Charge and matter densities are among the most funda- 
mental properties of atomic nuclei. Nuclear charge distri- 
butions have been intensively studied with elastic electron 
scattering and via the spectroscopy of X-rays from muonic 
atoms (see e.g. P0EED- These experiments profit from 
the full understanding of the electromagnetic interaction. 
Analyses of the distributions in the frameworks of differ- 
ent models have extracted characteristic parameters like 
charge radii, skin thicknesses, or the central depression 
of the charge density with high precision [^]|fj]. However, 
all these properties are only related to the distribution of 
the protons in the nucleus. The electromagnetic interac- 
tion provides only very limited information on the neu- 
tron. Therefore, the extraction of neutron distributions, 
respectively nuclear matter distributions (i.e. the sum of 
proton and neutron density) is much less straight forward. 
Some results for specific single neutron orbits have been 
obtained with elastic magnetic electron scattering, mak- 
ing use of the magnetic form factor of the neutron \7\ 
However, most experimental methods use hadron in- 
duced reactions such as for example a, proton, pion or 
kaon scattering from nuclei. The analysis of such reac- 
tions requires scattering theories for strongly interacting 
particles, with all their uncertainties. An overview over 
the different methods can be found in 151. The situation 
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is such, that a systematic analysis of the nuclear matter 
distributions is still missing and in many cases the spread 
between the results obtained with different probes is still 
larger than the predicted differences in the proton and 
neutron distributions, which are on the order of 0.05 fm - 
0.2 fm for the rms radii of heavy nuclei |10| . 

The present paper summarizes the analysis of recent 
experimental results for the coherent photoproduction of 
7T° mesons in view of nuclear matter distributions. 
This reaction is particularly attractive as a complemen- 
tary method for the study of nuclear matter distributions 
of stable nuclei. As discussed below, in the energy region of 
interest, protons and neutrons contribute identically with 
the same amplitude. Furthermore, in contrast to hadron 
induced reactions it is not restricted to the nuclear surface 
but probes the entire nuclear volume. In this sense it is the 
ideal reaction to test the matter distribution in the bulk of 
a nucleus. A first attempt to determine nuclear mass radii 
with this method was made by Schrack, Leiss and Penner 
in 1962 ^21- However, at that time the achievable experi- 
mental precision was very limited and much inferior to the 
results from hadron induced reactions which profit from 
the large cross sections characteristic for the strong inter- 
action. Subsequently, an attempt was made to measure 
nuclear matter radii via the coherent photoproduction of 
p° mesons |13j . However, due to the experimental diffi- 
culties and the previously not well developed theoretical 
tools for the extraction of the form factors from coherent 
photoproduction reactions, the method was never system- 
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atically explored. On the experimental side the progress 
made in accelerator and detector technology during the 
last fifteen years has considerably enhanced our possibil- 
ities for the study of photon induced meson production 
reactions. In particular, the new generation of quasi con- 
tinuous beam electron accelerators like CEBAF in New- 
port News, ELSA in Bonn, ESRF in Grenoble, MAMI in 
Mainz and SPring8 in Osaka are all equipped with state- 
of-the-art tagged photon facilities and highly efficient de- 
tector systems, most of them with almost 4tt solid an- 
gle coverage. Profiting from this developments, recently 
much more precise cross section measurements for coher- 
ent 7T° photoproduction from carbon, calcium, niobium 
and lead nuclei became available At the same time, 
progress was also made in the theoretical understanding 
of these reactions. Modern calculations treating the ele- 
mentary process in the framework of unitary isobar mod- 
els, incorporating final state interaction (FSI) in distorted 
wave impulse approximation, and including in-medium ef- 
fects of the Z\-isobar with phenomenological self-energies 
have become available 52 • The purpose of this paper is 
thus to explore the currently achievable accuracy in the 
determination of the nuclear mass distributions from co- 
herent 7T° photoproduction. This is also done in view of 
the possibility of further improvements in the data qual- 
ity. Although the data from [2], which are the basis of 
the present analysis, are the most precise results for co- 
herent 7r°-photoproduction from heavy nuclei available so 
far, they have been measured with an early stage of the 
TAPS detector [T311TC] covering only w 20% of 4tt. This 
resulted in typical detection efficiencies for ir° mesons on 
the order of only a few per cent and an imperfect sup- 
pression of incoherent contributions from excited nuclear 
states. With the availability of Att detector systems like 
the combined Crystal Ball/TAPS setup, operating now 
at the MAMI accelerator, more precise data will become 
available in the future. 

In this paper, we will limit the discussion of the mass 
distributions to the extraction of rms radii from the slope 
of the form factors and to the extraction of diffraction radii 
and surface extensions in the Helm model. However, the 
analysis of the full matter distributions in the same model 
independent way as for the charge distributions should 
become feasible for the next generation of experiments. 



2 Coherent 7r° photoproduction from nuclei 

Coherent photoproduction of 7r°-mesons from a nucleus 
with mass number A is the reaction 



7 + A(gs) A(gs) + ir° 



(1) 



where A{gs) is a nucleus in its ground state. It can be 
experimentally separated from breakup processes, where 
nucleons are removed from the nucleus, via its character- 
istic two-body kinematics fTJ]. The theoretical treatment 
is much more straight forward than for incoherent pion 
production reactions, since for initial and final state only 
ground state properties of the nucleus are needed. 



In general, the isospin structure of the elementary pro- 
cess of 7T° photoproduction from the nucleon is given by 

A{ 1 p-,^p) = +^-A V3 + ^-{A IV -A IS ) (2) 

A{ in - 7i° n?) = +y| A V3 + ^-{A IV + A IS ) , 

where A IS , A IV , and A V3 are the isoscalar, isovector, and 
total isospin changing parts of the total amplitude. How- 
ever, at incident photon energies in the range of interest 
in this work (200 - 350 MeV) the reaction is completely 
dominated by the photo excitation of the 4(1232) reso- 
nance (see e.g. [17]). Since this is an isospin 1 = 3/2 state 
only the isospin changing vector component A V3 can con- 
tribute, so that for the A excitation 



A(-yp — » tt°p) = A(jn — » TT°n) 



(3) 



Detailed investigations of coherent and breakup photopro- 
duction of 7T° mesons from the deuteron [2T)ll21lll7| have 
confirmed that the elementary cross sections for protons 
and neutrons are equal. This means that apart from small 
background contributions (nucleon Born terms) protons 
and neutrons contribute with the same amplitude to co- 
herent 7T° photoproduction from nuclei, so that this re- 
action is indeed sensitive to the distribution of nucleons 
rather than to the distribution of charge in the nucleus. 

In the most simple plane wave impulse approximation 
(PWIA) the coherent cross section from spin zero nuclei 
can be written as |14llll| 



da pwia 
dQ 



(E 1 , Q„) 



N 



dQ 



f(£*,e*)i^) (4) 



ddNS 



dQ 



±(E%,<%) = -g|.F 2 (i^)|W(0;), (5) 



where E~ i and n are incident photon energy and pion 
polar angle in the photon-nucleus cm-system, mjv is the 
nucleon mass, q(E~,, 6^) the momentum transfer to the nu- 
cleus, and F(q) the nuclear mass form factor. The total 
energy \fs of the photon-nucleon pair, the photon energy 
and momentum E*, k*, and the pion angle and momen- 
tum 0* , q* in the photon-nucleon cm-system can be evalu- 
ated from the average momentum pjv of the nucleon in the 
factorization approximation = q{A — 1)/2A. The spin- 
independent elementary cross section da^s/dQ is calcu- 
lated from the isospin average (for I nuclei weighted 
with N,Z) of the standard Chew-Goldberger-Low-Nambu 
(CGLN) amplitude Ti f^H] taken from ^H]. The extraction 
of the form factor F(q) from the differential cross sections 
in this approximation is straight forward and used below 
for a first approximative determination of the mass radii. 

It is well known, that final state interaction effects can 
have a significant impact on the pion production cross 
sections. The available model calculations have been com- 
pared to the measured differential cross sections in 
The distorted wave impulse approximation in momentum 
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PWIA (Kamalov) 
DWIA (Kamalov) 



Fig. 1. Total cross section for coherent n° photoproduction 
from carbon, calcium and lead (11) compared to the model 
results from |14|. 



space with additional A in-medium effects by Drcchscl et 
al. 25 gave the best agreement with the data. The results 
from this model are therefore used in the present work for 
the correction of FSI effects in the extraction of the nu- 
clear form factor. It should be emphasized that the model 
[H) was not adjusted to the nuclear data under discus- 
sion. The free parameters for the Z\-nucleus phenomeno- 
logical self energy were fitted to coherent ir° photoproduc- 
tion from 4 He and not modified for the heavy nuclei. The 
typical agreement between data and model results for the 
total cross sections is shown in fig. ^ It is quite good for 
carbon and calcium, but less so for lead. However, also 
for lead important features for this analysis, like the po- 
sition of the diffraction minima, are very well reproduced 
(see 22 for a detailed discussion). For the present analy- 
sis of form factors only the relative shifts of the position 
of the minima and cross section ratios between the PWIA 
calculation and the full model from 2H are used. Such cal- 
culations are presently not available for the nucleus 93 Nb. 
However, the corrections for the shift of the minima and 
the slope at small q can be approximated from lead by 
normalizing the position of the first diffraction minimum. 
Strictly speaking, also the above PWIA approximation for 
spin zero nuclei is not valid for this J — 9/2 odd-even nu- 
cleus. However, since it was shown in jllj that the mea- 
sured cross sections scale in the same way as for the spin 
zero nuclei (the contribution from the odd nucleon is not 
significant) we have kept it in the analysis. The systematic 
uncertainty is of course larger than for the other nuclei. 

An example for the influence of FSI on the angular 
distribution from lead 1 for an incident photon energy of 
290 MeV is shown in fig. H where the PWIA and DWIA 
calculations for the differential cross section are plotted 
versus the momentum transfer to the nucleus. The main 
effects at these incident photon energies are a reduction 
of the plane wave cross section and a shift of the position 
of the diffraction minima due to the pion - nucleus poten- 
tial. The latter is more pronounced for the higher order 
maxima. 



1 Note: differential cross sections da/dO are given as function 
of q throughout this paper, which for fixed i? 7 is a unique 
function of 8 n . This has the advantage that the position of the 
diffraction minima is approximately independent on E 1 . 



1 1.5 

q[fm" 1 ] 

Fig. 2. PWIA and DWIA calculations for coherent 7r° photo- 
production from lead for an incident photon energy of _E 7 =290 
MeV I14lll| . The shift of the diffraction minima due to FSI is 
indicated. 



For the interpretation of the form factor results it is 
of interest if final state interaction effects, in particular 
pion absorption, are so strong that like in hadron induced 
reactions effectively only the nuclear surface is tested. 

An indication for the strength of these effects can be 
obtained from the scaling of the cross sections with mass 
number in a similar way as in |22U23| for incoherent me- 
son production (n, rf) from nuclei. The production cross 
sections without FSI for incoherent processes scale with 
the number of nucleons A. The measured cross sections 
scale like A 2 / 3 , which implies strong FSI corresponding 
to a small mean free path of the mesons so that effec- 
tively only the nuclear surface contributes. The cross sec- 
tion for coherent ir° production in plane wave without FSI 
scales like A 2 , since the amplitude is proportional to A. If 
only the surface would contribute one would expect a scal- 
ing with A 4 ^ 3 . The observed scaling for the total coherent 
cross section does not even reach A 2 / 3 [23- However, one 
must take into account the influence of the sin 2 (0)F 2 (q) 
term in the PWIA cross section, which contributes to the 
A-dependence of the total cross section (see eq. ©). The 
nuclear form factor at q normalized to the position of the 
first diffraction minimum, q = qi(A), is almost indepen- 
dent on A in the region before the first diffraction mini- 
mum (see fig. El Consequently, the mass number scaling of 
the coherent cross sections can be obtained by fitting the 
ratio Rpwia (see eq. (J7J)) at a constant value of q/qi < 1 
with the ansatz 



Rpwia(A) cx A c 



(0) 



4 



B. Krusche: Nuclear mass form factors 



Q_ 

DC 

10 



-1 




atomic mass number A 




1.5 



1 

200 225 250 275 300 325 

E y [MeV] 

Fig. 3. Upper part: ^-dependence of Rpwia at q = 0.5gi for 
incident photon energies of 210, 230, 255, 280, and 305 MeV 
(from top to bottom). Lower part: fitted coefficients a of the 
mass dependence. 



The result for q = qi(A)/2 (i.e. approximately in the 0-th 
maximum 2 of the differential cross section) is shown in 
fig. At the lowest investigated incident photon energies 
around 200 MeV the scaling is very close to A 2 , indicat- 
ing almost negligible pion absorption. At higher incident 
photon energies FSI effects become more important. 



3 Extraction of mass radii 

3.1 rms-radii 

The determination of the root-mean-square (rms) radii 
requires the extraction of the nuclear form factor F(q 2 ) 
from the angular distributions. This is done in three steps 
as demonstrated in fig.r^for 40 Ca for one range of incident 
photon energy. In the first step the form factor is extracted 



2 Note: the maxima of the differential cross section are la- 
beled 0,1,2,... ,the maxima of the form factor 1,2,3,... . In both 
cases maximum 1 follows the first diffraction minimum. 
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Ca E =200 - 245 MeV 
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Fig. 4. Extraction of the nuclear form factor for 40 Ca. Trian- 
gles: measured differential cross section (axis labeling at right 
hand side). Open circles: PWIA approximation of form fac- 
tor from eqs. I0J,@. Filled circles: PWIA approximation after 
correction of angular resolution effects. Filled squares: after 
correction of FSI effects (DWIA approximation). 



from the data in plane wave approximation, i.e. all FSI 
effects are neglected. This is done in the same way as in 



[TT] 3 using 



F\q)\ 



PWIA 



Rpwia 



' N 



df2* 



with daNs/dd* from eq. J5J. 

In the next step a correction is applied for the finite 
angular resolution of the experiment. It was determined 
with Monte Carlo simulations that the detector response 
for 7T° mesons in the kinematical regime of interest corre- 
sponds to a Gaussian with a FWHM of 4° for the pion cm 
angles. The DWIA model calculation was folded with this 
response function and the data were corrected with the 
ratio of folded and unfolded calculation. This correction 
is only significant for niobium and lead. 

In the final step the ratio of PWIA and DWIA cross 
sections obtained from the model calculations is used to 
correct the form factor for FSI effects. The q dependence 



3 Note: eq. [7| is applied to the data in the finest possible 
binning of incident photon energies, defined by the resolution of 
the tagging spectrometer (roughly 2 MeV). The corresponding 
form factor results are then averaged over larger energy regions 
for statistical reasons. 
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(not the absolute values) of the P WIA and DWIA calcula- 
tions is in most cases similar for small g, nevertheless this 
correction introduces an additional model dependence into 
the analysis. We have therefore extracted the mass radii 
from the form factors with and without this correction in 
order to get an estimate for the systematic uncertainties. 

Once the form factor has been determined, the rms- 
radius can be extracted without further model assump- 
tions from the slope of the form factor for q 2 — > via 



F(q 2 ) = l-Y* rms + 0(q% (8) 
which is done in the usual way by fitting a polynomial 



JV 



i? (9 2 ) = E(- 1 )" a «'? 2n 



(9) 



n=0 



to the data. The rms-radius is then given by 

r rms = A/6ai/a , (10) 

where for a correctly determined form factor ao should 
be unity. This is of course not true for the form factors 
extracted in PWIA approximation without correction for 
FSI effects. However, also the form factors extracted in 
DWIA approximation differ in most cases somewhat from 
F(q 2 ) — ► 1 for q 2 — ► (see figs. 16170 . This can be due 
to systematic uncertainties in the absolute normalization 
of the measured cross sections (up to 10% overall, up to 
3% relative between different nuclei, see ref. |U]) or due 
to an imperfect correction of the FSI effects (systematic 
uncertainties in the models). Therefore ao is kept free in 
all fits to account for these effects. 

A final remark must be made to the comparison of the 
mass form factors, extracted in this way, to the nuclear 
charge form factors. In the latter case the charge distribu- 
tion of the nuclei is tested. Due to the charge form factor 
of the proton the distribution of (point-like) protons F pc 
in the nucleus is given by 



F pc {q) 



F ch {q) 
FWqj 



(11) 



where F ch is the nuclear charge form factor and F p h is 
the charge form factor of the proton for which we take the 
dipole form factor 

2 \ - 2 

q 



jpch 



1 



18.234fm" 



(12) 



On the other hand, coherent pion production, which pro- 
ceeds through the excitation of the nucleon to the A res- 
onance, is testing the distribution of point-like nucleons 
in the nucleus. Therefore for the comparison of the mass 
and charge radii the rms charge radius r^ ms = 0.862 fm of 
the proton was subtracted in quadrature from the nuclear 
charge rms-radii r c ^ ms to give the rms-radius r^ na for the 
distribution of point-like protons 



rms 



ch )2 
rms J 



( r P \2 
\i rrns J 



(13) 



For the comparison of skin thicknesses the charge form 
factors were divided by the proton charge form factor. 



3.2 Form factors in the Helm model 

The extraction of the rms radii from the slope of the form 
factor has the advantage that no model of the form fac- 
tor itself is needed. However, models of the form factor 
which relate for example the radii to the position of the 
diffraction minima allow a much better control of system- 
atic effects like the DWIA corrections. Furthermore, addi- 
tional information can be gained from such models. The 
rms radius alone has for example no information about 
the extension of the surface zone of the nuclei. A good 
example are the charge rms radii of 40 Ca and 48 Ca, which 
are almost identical. However, the actual charge distribu- 
tions are by no means identical. The nucleus Ca has a 
larger core region of almost constant density but a smaller 
surface region, where the density drops from 90% to 10 5 . 
For the model dependent analysis it is convenient to use 
Helm's model, which is known from the analysis of elec- 
tron scattering data (see |3|H]) and allows to extract nu- 
clear extension parameters in a transparent way. In this 
model the nuclear density is parameterized [S] by the con- 
volution of a hard sphere distribution with a Gaussian. 
The form factor Fjj is then simply given by the product 
of the form factor of the hard sphere, Fh s , with that of 
the Gaussian, which, again is a Gaussian, Fq. 



F, 



H 



Fg ■ Fh 



F G = exp(-{q<j) 2 /2) 



(14) 
(15) 



Fhs — 



sin{qR d ) \ 3 . 

cos(qR d ) = ——j 1 {qR d ) 

qttd J qtid 

Here, q is the momentum transfer and j\ the first order 
spherical Bessel function. R d is the so-called 'diffraction 
minimum sharp radius' (dms-radius) and the width of the 
Gaussian a is approximately related to the 10%-to-90% 
surface width tn via 



t H = 2.54cr 



(17) 



In this model, the zeroes of the form factor are deter- 
mined by the zeroes of the Bessel function, which implies 
a straight forward relation between the dms radius and 
the momentum transfers q^ i=l,2,3,... in the i-th mini- 
mum of the form factor. The rms radius, which integrates 
over the core region and the surface zone, is then related 
to the dms radius and the a of the distribution by 



(18) 



The a of the Gaussian can be extracted e.g. from the po- 
sition q m of the first maximum (i.e. the maximum after 
the first minimum) of the form factor and its magnitude 
F(q m ) via |S] 




2 , 3ji(q m R d ) 
—prln- 



m J 



(19) 



The comparison of the dms mass radii to the dms 
charge radii from electron scattering requires no correc- 
tion for the proton charge radius since multiplication of 
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the form factor eq. (|14|) with the proton dipole form fac- 
tor has no influence on the position of the zeroes, which 
are still determined by the Bessel function. This is not 
the case for the width of the Gaussian which is expected 
to be larger in case of the charge form factor due to the 
contribution of the proton charge form factor. 

4 Results 

In the following we will discuss the form factors extracted 
for 12 C, 40 Ca, 93 Nb and ™ at Pb 4 from the pion production 
data. We will first discuss the overall features of the mass 
form factors (without DWIA corrections). In the second 
subsection we discuss the extraction of root-mean-square 
radii from the slope of the form factors and in the final 
subsection the interpretation of the form factors in the 
framework of the Helm model. 



4.1 Form factors in PWIA 

The form factors have been extracted in PWIA approx- 
imation from the cross section data for incident photon 
energies from 240 - 300 MeV with eqs. gj , (0 , Q . At 
this stage only the overall normalization was corrected 
for DWIA effects. For this purpose the form factors were 
fitted at low momenta with the Taylor series eq. © and 
the data where divided by the fitted a Q coefficient. The 
results are summarized in fig. [SJ where the form factors 
are plotted versus qRd which is calculated as: 

qR d {A)=q-X 1 /q 1 (A) (20) 

where q\ (A) is the momentum transfer in the first diffrac- 
tion minimum and Xi = 4.495 is the first zero of the 
Bessel function. In this plot the positions of all diffraction 
minima fall on top, as one would expect as long as the 
Helm model is valid, and even more the form factors for 
Ca, Nb and Pb are almost identical and agree up to the 
first diffraction minimum with the form factor of a hard 
sphere. Only the carbon form factor shows the effect of a 
finite surface region in the range before the first minimum. 
The positions of the higher order minima and maxima do, 
however, not agree with the form factor of a hard sphere 
(grey curve), since the shift caused by the FSI effects is 
not the same for the different minima (in case of lead 
roughly 3% for the first minimum and 20% for the 5-th 
minimum). However, if the zeroes of the hard sphere form 
factor are adjusted to these shifts (full curve) the data are 
quite well described by the hard sphere form factor. This 
is surprising since due to the effects from the finite surface 
region the Helm model predicts a significant faster fall-off 
of the form factors in the region beyond the first diffrac- 
tion minimum. To demonstrate this effect the form factors 

4 Note: In all cases targets with natural isotope composition 
have been used in the experiments However, natural nio- 
bium is mono-isotopic, and in case of carbon and calcium the 
admixture of other isotopes than 12 C and 40 Ca is negligible at 
the current level of precision. 



of a hard sphere convoluted with a Gaussian with a— 0.75 
fin and the R4 radii of Ca (dotted) and C (dashed) are 
also shown in the figure. The predicted magnitude of the 
form factor in the first maximum decreases strongly with 
decreasing mass number, but the experimental results are 
practically identical for all nuclei. It will be discussed in 
subsection 14.31 that the effect of the surface thickness is 
canceled to a large extent by DWIA effects, so that the 
extraction of surface parameters requires a careful correc- 
tion of the pion distortion effects. 




Fig. 5. Mass form factors extracted in PWIA approximation 
(see text). Grey curve: hard sphere form factor (eq. 1161 h solid 
curve: hard sphere form factor with zeroes adjusted for DWIA 
effects, dotted: hard sphere (Rd of Ca) and Gaussian with 
0= 0.75 fm, dashed: hard sphere (Rd of C) and Gaussian with 
cr=0.75 fm. Insert: region before first minimum in linear scale. 



4.2 rms-radii 

For the form-factor model independent extraction of the 
rms-radii the results with and without correction of FSI ef- 
fects (DWIA, PWIA approximation, see fig. El were used. 
The DWIA approximation is of course expected to give 
more realistic results. However, also in this approxima- 
tion the form factors do not exactly approach unity for 
<Z^0. 

The size of these effects is summarized in tabled where 
the overall normalization constant a (see eq. (0) is listed 
for the PWIA and DWIA approximations for two different 
regions of incident photon energy. In all cases averages for 
fits with polynomials of degree iV=2,4 are given, however 
the iV=2,4 results differ only by a few per cent. The DWIA 
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E y =200 - 245 MeV 




0.2 0.4 0.6 0.8 1 
q 2 [fm" 2 ] 

Fig. 6. Fit of the form factors in PWIA and DWIA with poly- 
nomials of degree N=4 (solid lines) and N=2 (dashed lines) 
for incident photon energies between 200 MeV and 245 MeV. 



approximation brings the overall normalization closer to 
unity, in particular for the higher incident photon energies. 

Since the rms-radii depend only on the slope of the 
form factor but not on the absolute normalization (as long 
as a is kept a free parameter, see eqs. |JSJ-|QDJ), it is a pri- 
ori not clear how large the influence of the FSI corrections 
is. In order to get some estimate, results from the PWIA 
and DWIA approximation are compared. Typical fits of 
the data are shown in figs. 16171 The results of the fits for 
the rms-radii calculated from eq. IjlUI) are summarized in 
fig. They are plotted as function of the upper limit of 
the fit range, which allows to judge the stability of the fits. 
Fits have been done with polynomials of degree N=2,4. 
The fits with N—A become only stable when rather large 
ranges of q are fitted, however in that limit they approach 
the Af=2 results. The results from PWIA and DWIA agree 
in most cases better than within 10%. However, the radii 



E y =245 - 290 MeV 




0.2 0.4 0.6 0.8 1 
q 2 [fm" 2 ] 

Fig. 7. same as fig. |S| for photon energies between 245 MeV 
and 290 MeV. 



extracted in DWIA approximation are in all cases system- 
atically smaller than those from the PWIA analysis, this 
effect increases with nuclear mass. This is expected, since 
the FSI effects tend to increase the slope of the form factor 
(they shift the first minimum to smaller q values, see fig. 



Table 1. Normalization constants of the form factor data. 
Listed is the a coefficient of the fitted Taylor series (see eq. 
©). The values are the average for fits with polynomials of 
degree iV=2,4. 





200-245 MeV 


245-290 MeV 


nucleus 


PWIA 


DWIA 


PWIA 


DWIA 




0.95 


1.20 


0.70 


1.10 


40 Ca 


0.83 


1.11 


0.55 


0.95 


93 Nb 


0.75 


1.07 


0.45 


0.92 


na *Pb 


0.72 


0.99 


0.42 


0.83 
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o PWIA, n=2 □ PWIA, n=4 
• DWIA, n=2 ■ DWIA, n=4 



a PWIA adopted 
a DWIA adopted 




Fig. 8. Fit results for the rms-radii extracted from eqs. JHJ- 
J 1UP . They are plotted as function of the upper limit of the 
fitted g-range for the PWIA and DWIA analyses and for fits 
with polynomials of degree N = 2, 4. They have been averaged 
over the two ranges of incident photon energies. The triangles 
represent the adopted values, the stars the corresponding rms 
charge radii from ref. £Q. 

|2J), which leads to an overestimate of the radii in PWIA 
approximation. The adopted values of the radii are com- 
pared in fig. [H] and in table to the nuclear charge radii 
extracted from electron scattering experiments. The un- 
certainties have been estimated from the variation of the 
results with the fit range and the agreement between the 
N=2 and N=4 results. 

Table 2. Mass rms radii (r rms ) extracted in PWIA and DWIA 
compared to charge rms radii (rfms) DP- The charge radii have 
been corrected for the proton charge form factor by subtracting 
in quadrature the proton rms radius (r p =0.862 fm). 



nucl. 




[fm] 


' rms 


r lr pc 

1 rms I I rms 




PWIA 


DWIA 


[fm] 


PWIA 


DWIA 


12 C 


2.28±0.10 


2.26±0.10 


2.32 


0.98 


0.97 


40 Ca 


3.22±0.10 


3.15±0.10 


3.35 


0.96 


0.94 


93 Nb 


3.96±0.10 


3.76±0.10 


4.22 


0.94 


0.89 




5.20±0.15 


4.90±0.15 


5.42 


0.96 


0.90 



In all cases the mass radii are somewhat smaller than 
the corresponding charge radii and this effect becomes 
larger for increasing mass. The effect lies between 2-6 
% for the PWIA approximation and between 3 - 10 % 
for the DWIA results (note that since the DWIA correc- 
tion for 93 Nb was approximated from Pb it has a larger 
systematic uncertainty). This result is in particular un- 



expected for lead, where nuclear models (see e.g. ^01) in- 
general predict somewhat larger rms radii for the neutron 
distribution than for the protons, so that one would to 
the contrary expect slightly larger mass rms radii. Possi- 
ble explanations include still uncontrolled effects from the 
correction of the FSI in DWIA and small contaminations 
of the coherent cross section with incoherent reaction pro- 
cesses (see sec. EJ. 



4.3 Results from the Helm model 

The most sensitive analysis of the form factor in the frame- 
work of the Helm model would be a fit of the data over the 
full measured q range. However, it is not straight forward 
to apply the FSI corrections, which include shifts of the 
position of the minima and modifications (reductions for 
most incident photon energies) of the magnitude, consis- 
tently over the full q range. Therefore, the direct deter- 
mination of Rd from the position of the minima and the 
determination of a from eq. Ijl9|l has the advantage that 
only specific features (shift of minima and modification 
of first maximum) of the DWIA corrections of the cross 
sections must be known. At this stage we do not attempt 
to fit the form factors beyond the first minimum, where 
the DWIA effects can be corrected in the way discussed in 
sec. 13.11 However, more refined treatments of the DWIA 
corrections are certainly possible and will be worthwhile 
when more precise data become available. Here the anal- 
ysis is done in two different ways. 

Table 3. Diffraction mass radii (Rd) and width of the Gaus- 
sian (a) determined from the fit of eq. 1141 to the data. The 
parameters are averages over the two regions of incident pho- 
ton energies (200 -245, 245 - 290 MeV). The uncertainties have 
been estimated from fit errors and the agreement between the 
two energy ranges. For comparison the parameters from charge 
distributions R d , a d are also listed. They represent the aver- 
ages of the results given in pj (since 93 Nb is not analyzed in 
instead the average for the isotopes 92,94 Zr is given in brack- 
ets). 



nucleus 




a[fm] 


R d [im] 


<7c [fm] 


12 C 


2.30±0.10 


0.7 - 1.2 


2.44 


0.8 


40 Ca 


3.65±0.30 


0.6 - 1.2 


3.79±0.04 


0.92 


93 Nb 


5.10±0.10 


0. - 1. 


(5.04±0.04) 


(0.98) 




6.65±0.05 


0. - 1. 


6.66±0.04 


0.93 



In the first analysis the form factors in the region be- 
fore the first diffraction minimum extracted in the DWIA 
approximation as discussed in sec. 13.11 have been fitted 
with eq. Ijl4|) multiplied with an overall normalization con- 
stant. The fits are shown in fig. E| for two ranges of inci- 
dent photon energies. The results for the fit parameters 
are summarized in tab. For the dms radii quite good 
agreement is found with the charge radii extracted from 
electron scattering experiments. Since the fitted g-range 
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CNJ 



10 




10 



E =245 - 290 MeV 




q 2 [fm" 

Fig. 9. Fit of the Helm model (eq. dJ) to the form factors 
in DWIA approximation for two different ranges of incident 
photon energies. 



is small for the heavy nuclei, the fits are not sensitive to 
the width of the Gaussian therefore in case of Nb and Pb 
values between zero and unity result in similar fit qualities. 

A more precise determination of the radii uses the po- 
sition of all observed diffraction minima. The results for 
Ca and Pb are shown in figs. HOlfD For the first minimum 
of Ca and the first and second minima of Pb the positions 
were individually determined in bins of 20 MeV incident 
photon energy. Since no systematic trends where observed 
only the average over all bins is shown in the figures. Due 
to the limited statistics the higher minima where only de- 
termined for the full energy range from 245 - 365 MeV 
without considering possible energy dependent shifts. Al- 
though we have assigned rather conservative uncertainties 
to the positions of the higher minima, their large lever arm 
leads to fairly precise values of the corresponding values 
of the radii. 



Ca(Y,7t°)Ca 



E = 245 - 365 MeV 



4.5 



E 

^4 



It 



3.5 



'A | a , 

+ Ta 



12 3 4 

number of minimum 



♦ +4 



0.5 



1 



1.5 



2.5 



momentum transfer q[fm" ] 



Fig. 10. Main plot: position of the diffraction minima for 
40 Ca. Insert: Extracted dms radii. Open squares: raw values, 
filled circles: after correction for FSI, triangles: charge dms- 
radii from electron scattering 0- 



Pb(y,7i°)Pb 



E = 245 - 365 MeV 



10 



^ 10^ 



10 



E 



* A * ^ u 



1 2 3 4 5 6 

number of minimum 



!»,((( 



lit 



li 



0.5 1 1.5 2 2.5 

momentum transfer q[frrf 1 ] 

Fig. 11. Same as fig. ED for Pb. 



The radii have been calculated from the position of the 
minima via (see eq. Hlfij l): 



R 



(i) 



X i /q l 



(21) 



where qi is the position of the i-th minimum, Xi the i- 

(i) 

th zero of the Bcsscl function and R d the corresponding 
value of the dms radius. The results are shown in the in- 
serts of the figures and summarized in fig. 1121 

The values extracted from the raw data with eq. (|21J) 
(open squares in the figures) show a systematic trend as 
function of the number of the minima. However, this trend 
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is completely eliminated and almost perfect agreement 
with the corresponding charge radii is obtained after cor- 
rection of the FSI effects (filled circles in the figures). The 
shift of the position of the minima due to FSI effects was 
determined from a comparison of the position of the min- 
ima in the PWIA and DWIA calculations as discussed in 
sec. (see fig-EJ- The good agreement between the values 
extracted from the different minima after correction nicely 
demonstrates that the FSI effects are well under control 
at least as far as the position of the minima is concerned. 




D 1 2 3 4 5 6 

number of minimum 

Fig. 12. Extracted values of the dms radii. Open squares: raw 
values, filled circles: after correction for FSI, triangles: charge 
dms- radii from electron scattering ,5.. In case of Nb the open 
triangles correspond to the average of the charge dms radii for 
92 ' 94 Zr (values for 93 Nb are not available). The full and dashed 
lines indicate average and uncertainty of the radius for each 
nucleus. 



For the charge distributions a systematic deviation of 
the dms radii extracted from the first and second diffrac- 
tion minium has been found |2K| , which was related to the 
depression of the central density of the charge distribution 
in nuclei. It would be therefore very interesting, to investi- 
gate if such a central density suppression survives also into 
the mass density. However, the precision achieved so far, 
does not yet allow to answer this question. The extracted 

(R ( d 1) - R { P) differences are 0.10 ± 0.12 fm, 0.16 ± 0.13 
fm, and 0.00±0.15 fm for 40 Ca, 93 Nb and Pb. This means 
that they are consistent with zero but also with the small 
differences observed for the charge radii (on the order of 
0.1 fm for lead). 

Finally the first maximum of the form factors was used 
for a more precise extraction of the width of the Gaussian 
with eq. id. This requires not only the determination 
of the position q m of the maximum but also the absolute 
value of the form factor F(q m ), which is strongly influ- 



enced by the FSI effects. In order to minimize the model 
dependency, the correction was done in the following way. 
As in sec. 14.11 the form factors were fitted with Taylor 
series and the absolute normalization was obtained from 
the condition a =l. The results are shown in fig. El left 
hand side. Position q m and magnitude F(q m ) of the first 
maximum were determined from these data. The FSI cor- 
rection for the position was obtained in the same way as 
for the position of the minima. The correction of the mag- 
nitude was also obtained from a comparison of the differ- 
ential cross sections calculated in PWIA and DWIA. This 
is shown on the right hand side of fig. El In the picture, 
the position of the first minimum and the magnitude be- 
fore the first minimum are normalized for the DWIA cal- 



Table 4. Determination of a from the first maximum of the 
form factor (see eq. 1191 ). q m and F(q m ) are position and mag- 
nitude of the first maximum, q^ r and F cor (q m ) the same after 
correction for FSI effects, (q in [fm~l], a in [fm]). 



nucl. 


q m 


cor 


F(q m ) 


F cor (q m ) 


1m ltd 


a 


12 C 


1.95 


2.15 


0.079 


0.043 


5.18 


0.46 


40 Ca 


1.30 


1.41 


0.076 


0.058 


5.33 


0.53 


nai Pb 


0.78 


0.825 


0.0776 


0.071 


5.49 


0.67 



culation so that the cross sections agree in the q range 
before the first minimum. The additional FSI correction 
for the magnitude of the first form factor maximum is then 
obtained from the square root of the ratio of the two cal- 
culations at this q value. In this way, again only ratios of 
the model results enter into the correction factors. In the 
absence of a good DWIA calculation for Nb no analysis of 
the Gaussian width was attempted for this nucleus. The 
parameters of the first maximum and the deduced values 
for the Gaussian width a are summarized in tab. ^ 

It is obvious from figure and table fj] that the rela- 
tive FSI corrections (i.e. the difference of the FSI effects 
between the q range before the first minimum and in the 
first maximum) increase with decreasing mass number, 
they are almost negligible for lead but very significant for 
carbon (note that the absolute FSI effects are of course 
much larger for lead than for carbon). Their effect is, that 
for the form factors normalized to unity for q — > the 
magnitude of the first maximum is overestimated for the 
light nuclei in the PWIA analysis. And this effect cancels 
the influence of the finite Gaussian width on the magni- 
tude in the first maximum. 

The results for the dms radii and the Gaussian width a 
are summarized in tab.[S]and compared to the correspond- 
ing values for the charge distributions. The mass dms-radii 
are in excellent agreement with their charge counterparts, 
though the almost exact agreement for Ca and Pb is cer- 
tainly by chance, given the associated uncertainties. Thus 
also for the diffraction radii there is no indication that the 
mass radii would be larger than the charge radii. 
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1 2 

q[fm] 

Fig. 13. Left hand side: position and height of the first max- 
imum of the form factors for carbon, calcium and lead. Right 
hand side: PWIA (dashed curves) and DWIA (full curves) cal- 
culations for the differential cross sections for the same range 
of incident photon energies as taken into account in the data 
analysis (200 - 260 MeV for Ca and Pb, 255 - 330 MeV for C). 
Position of the first minimum and the absolute values before 
the first minimum of the DWIA calculation are normalized to 
the PWIA calculation. 



Table 5. Parameters of the Helm model from the positions of 
the diffraction minima and position and magnitude of the first 
maximum. The values for Rd are the average values from fig. 
1121 The results for a are from tab. 2] Rd an d o"c are the re- 
spective values for the nuclear charge distributions taken from 
|2] (the latter corrected for the finite extension of the proton 
charge distribution) . 



nucl. 


Rd[{m] 


cr[fm] 


| H5M 


<7c [fm] 


12 C 


2.41±0.09 


0.46 


2.44 


0.67 


40 Ca 


3.78±0.05 


0.53 


3.79±0.04 


0.80 


93 Nb 


5.09±0.05 




5.04±0.04 






6.66±0.07 


0.67 


6.66±0.04 


0.81 



For the comparison of the width a with the corre- 
sponding width <j c of the charge distributions one must 
take into account that the values for the charge distribu- 
tions include the effect of the finite extension of the proton, 
while the values extracted from pion photoproduction re- 
fer to point like nucleons. Therefore a c has been corrected 
for this effect with the dipole form factor of the proton. 
However, even after this correction, the values extracted 
for the mass distributions remain systematically smaller, 
in particular for the light nuclei (see tab.EJ. 



5 Summary and Conclusions 

Recent data for coherent photoproduction of it" mesons 
have been analyzed in view of nuclear mass distributions. 
The Helm model was used to extract diffraction-minimum 
sharp radii from the positions of the diffraction minima 
and skin thicknesses from the magnitude and position of 
the first maximum of the form factors. After corrections 
for FSI effects, the mass dms-radii are in excellent agree- 
ment with the dms-radii of the charge distributions of the 
nuclei. A stringent control of systematic effects (e.g. the 
FSI corrections) on the dms-radii is possible via the com- 
parison of the results for different order minima of the 
same nucleus. No systematic trends have been observed. 
The uncertainties of the dms radii (3.7 % for 12 C and « 
1 % for all other nuclei) have reached a similar level of 
precision as the corresponding charge radii. However, the 
achieved precision is just not sufficient to exclude or estab- 
lish a central depression of the mass distributions from the 
comparison of the positions of first and second diffraction 
minimum. 



Table 6. Comparison of the extracted rms mass radii with 
charge rms radii (all values are in [fm]). Rd'. mass dms radii (see 
table|2J, a: Gaussian width of the mass distributions (see table 
[HJl (value in brackets for Nb interpolated from other nuclei), 
rf^j : rms mass radii calculated with eq. (1181 from Rd and 
fj r rms, r ?ms'- rm s radii extracted from the slope of the form 
factors in PWIA and DWIA approximation (see table 0, r c rms 
charge rms radii with the proton charge radius subtracted in 
quadrature. 



nucl. 


Rd 


a 


Helm 
' rms 


r PW 
1 rms 


r DW 1 

' rms 


r c 

' rms 


12 C 


2.41 


0.46 


2.03 


2.28 


2.26 


2.32 


40 Ca 


3.78 


0.53 


3.07 


3.22 


3.15 


3.35 


93 Nb 


5.09 


(0.60) 


4.08 


3.96 


3.76 


4.22 


nat Pb 


6.66 


0.67 


5.29 


5.20 


4.90 


5.42 



Root-mean-square radii of the mass distributions have 
been extracted in two different ways. In the first, the slope 
of the form factors for q — > was fitted with polynomials. 
A very conservative estimate of the typical uncertainty 
due to FSI corrections of the shape of the form factors 
follows from the difference between the results obtained 
in PWIA and DWIA, which is more important the heav- 
ier the nucleus. In a second approach, they have been ex- 
tracted in the framework of the Helm model via eq. (f 1 H|) 
from the dms radii and surface thicknesses. The results are 
summarized in table HO The largest discrepancies between 
the different analysis are in the range of 10 %. 

An unexpected finding is, that all results for the rms 
mass radii are slightly smaller than the corresponding rms 
charge radii, even after the latter have been corrected for 
the proton charge radius. Since at least for nuclei like 
208 Pb models predict slightly larger rms radii for the neu- 
tron distributions than for the proton distributions one 
would have expected the opposite. At the same time the 
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skin thicknesses extracted from the Helm model are also 
systematically smaller than their charge counter parts. In 
a sense, these two effects are consistent since smaller skin 
thicknesses combined with identical diffraction radii will 
lead to smaller rms-radii. However, as yet it is not clear 
if these effects are real or if they can be explained by so 
far not understood systematic effects in the data or in 
the model calculations used for the DWIA corrections. A 
possible explanation could be a small remaining incoher- 
ent background component in the cross section data. The 
angular distributions of incoherent tt° photoproduction in- 
volving excited nuclear states peak at larger angles than 
the coherent reaction. Therefore, such background, which 
so far is only suppressed by kinematical cuts, would tend 
to decrease the slope of the form factor at small q, it would 
enhance the magnitude of the first diffraction maximum, 
but it would not change the position of the diffraction 
minima. In such a scenario, rms-radii and skin thicknesses 
would be underestimated, but the dms-radii would not be 
effected. Clearly, further improvements on the experimen- 
tal side and for the model calculations are necessary and 
possible. On the experimental side, large improvements in 
the statistical quality of the data are possible with the 
now available 47r electromagnetic calorimeters. Further- 
more, improvements in the energy and angular resolution 
for the detection of the ir° decay photons will allow an 
even more stringent suppression of incoherent pion pro- 
duction reactions than achieved in The background 
situation will furthermore very significantly improve since 
with the large solid angle coverage the detection of de- 
cay photons from excited nuclear states will not only al- 
low to veto such events much more efficiently than previ- 
ously, it will actually allow a detailed investigation of the 
incoherent processes (which are very interesting in their 
own right) so that any remaining background can be sub- 
tracted. Furthermore, a more systematic investigation of 
coherent pion photoproduction from many nuclei in the 
framework of models is clearly desirable. 
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